Semi-analytical approach to short-wavelength dispersion and modal properties of 

photonic crystal fibers 
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We consider photonic crystal fibers made from arbitrary base materials and 
derive a unified semi-analytical approach for the dispersion and modal prop- 
erties which applies to the short-wavelength regime. In particular we calculate 
the dispersion and the effective index and comparing to fully-vectorial plane 
wave simulations we find excellent agreement. We also calculate asymptotic 
results for the mode-field diameter and the ^-parameter and from the latter 
we predict that the fibers are endlessly single mode for a normalized air-hole 
diameter smaller than 0.42, independently of the base material. 
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Photonic crystal fibers (PCF) are dielectric optical 
fibers with an array of air holes running along the full 
length of the fiber. Typically, the fibers employ a sin- 
gle dielectric base material (with dielectric function Eb = 
nf) and historically silica has been the most common 
choice. 1 ^ Recently other base materials have been stud- 
ied including chalcogenide glass, lead silicate glass, tel- 
luride glass, bismuth glass, silver halide, teflon, and plas- 
tics/polymers. The fabricated fibers typically share the 
same overall geometry with the air holes arranged in a 
triangular lattice and the core defect being formed by the 
removal of a single air-hole. There has been a major theo- 
retical and numerical effort to understand the dispersion 
and modal properties of especially silica-based PCFs, but 
the scale-invariance of Maxwell's equations^ cannot be 
applied directly to generalize the results to other base 
materials. The reason is that PCFs made from different 
base materials do not relate to each other by a linear 
scaling of the dielectric function, e(r) s 2 e(r), with s 
being a positive scalar. The increased focus on use of 
new base materials thus calls for a theory of PCFs with 
an arbitrary base material. 

This Letter offers a unified approach to the disper- 
sion properties which utilizes that the base material typ- 
ically has a dielectric function exceeding that of air sig- 
nificantly, Eb 3> 1. The short- wavelength regime is char- 
acterized by having the majority of the electrical field re- 
siding in the high-index base material while the fraction 
of electrical field in the air holes is vanishing, see e.g. 
Refs. 0. The calculation starts from the fully-vectorial 
wave-equation for the electrical field^ 



V x V x E(r) = e(r)— E(r). 

c z 



(1) 



For a fiber geometry with z along the fiber axis we have 
e(r) = e(x, y) and we look for solutions of the plane-wave 
form e*^ 2- "') with the goal of calculating the dispersion 
relation uj([3). The above discussion for Sb 1 suggests 
that we can approximate the problem by imposing the 



boundary condition that E is zero at the interfaces to 
the air holes. Since the displacement field D = eE is 
divergence free we have = eV • E + E • Ve w eV • E. In 
the bulk matrix material the latter equality is exact and 
the wave equation reduces to 



-V 2 E(r)=e fc — E(r), 



(2) 



which we solve with the boundary condition that E is 
zero at the interfaces to the air holes. For the boundary 
condition to be meaningful it is crucial that Sb S> 1 so 
that the fraction of electrical field in the air holes is van- 
ishing in the short-wavelength regime. The wave prob- 
lem has now become scalar, but it should be emphasized 
that the approach is very different from the usual scalar 
treatment^ and recent applications to PCFs in the short- 
wavelength regim o 2 i 4 i 6 which take the electrical field in 
the air holes into account. Obviously, the scalar problem 
posed by Eq. J5J) is separable and formally we have that 



to 



-a 2 



(3) 



where fl z = c[3/nb is the frequency associated with the 
longitudinal plane-wave propagation, fl xy = 7 x chr x jnb 
is the frequency associated with the transverse confine- 
ment/localization, and 7 is a corresponding dimension- 
less and purely geometrical number, which only depends 
on the normalized air-hole diameter d/A. From Eq. (0) 
it follows that 7 is an eigenvalue governed by a scalar 
two-dimensional Schrodinger-like equation 



-A 2 {d 2 x + d 2 M{x,y) =j 2 i>(x,y), 



(4) 



with the scalar function ip being subject to hard-wall 
boundary conditions at the interfaces to the air-holes, 
i.e. ip = in the air holes. 

The developments in computational physics and en- 
gineering have turned numerical solutions of partial dif- 
ferential equations in the direction of a standard task. 
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Fig. 1. Panel (A) shows the geometrical eigenvalues j 2 
for the fundamental core (c, 1), the second order core 
(c, 2) and the fundamental cladding (cl) modes versus 
normalized air-hole diameter d/A, Panel (B) shows the 
corresponding ^-parameter, and Panel (C) the effective 
mode-field radius of the fundamental core mode. The 
data-points are obtained from finite-element simulations 7 
of Eq. J3J and the dashed lines are guides to the eyes. The 
gray region indicates the endlessly single-mode regime 
with Vpcf < t- 



Here, we employ a finite-element approach^ to numeri- 
cally solve Eq. Q and calculate 7 2 versus d/A. Panel (A) 
in Fig.^summarizes the results for the fundamental core 
mode (see inset in Fig. |2J), the second-order core mode, 
and the fundamental cladding mode (see inset in Fig. - 
For the core modes the problem has been truncated 
by considering a sufficiently large domain with Dirich- 
let boundary conditions while for the cladding modes we 
have considered the unit-cell with periodic boundary con- 
ditions which for symmetry reasons can be formulated in 
terms of Neumann conditions^ As the normalized air- 
hole diameter d/A is increased the localization becomes 
more tight and as expected the eigenvalue increases. For 
the fundamental core mode the dashed line shows a third- 



7c 2 ~ a/ + SB - + 



(5) 



with sf = 2.67, SB = 12.51, ^ = -9.45, and 9 = 13i 
being fitting parameters. 

The recently proposed ^-parameter— Vpcf = A(/3 2 
fici) 1 ^ 2 becomes 



lim Vpcf = ' ' ~ 1 



7d - % 



(6) 



and the numerical results shown in Panel (B) agrees 
nicely with the short-wavelength asymptotic limit of re- 
cent simulations on silica-based PCFs. 9 The endlessly 
single-mode regime, 2 defined by Vpcf < ^r 1 exists for 
d/A < 0.42 independently of the base material. For 
more detail on the modal cut-off see Ref. @ and refer- 
ences therein. Panel (C) shows results for the effective 
mode-field radius^ w e s calculated from A c s = irw 2 s with 
the effective area given by 

Jdxdy \t/j(x,y)\ Jdx'dy' \ipix' ,y')\ 2 
km A cff = '— -± '-. (7) 

A<A / dxdy \ip(x,y)\ 
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Fig. 2. Effective index n e ff of the fundamental core mode 
versus normalized wavelength A/ A for holey fibers with 
normalized air-hole diameter d/A = 0.4 and varying 
base material. The dashed lines are the predictions of 
Eq. JSJ and the data points are results of fully-vectorial 
plane-wave simulationsii The inset shows the fiber ge- 
ometry and the fundamental core eigenfunction ip c with 
7 2 = 7.0506, obtained with the aid of a finite-element 
simulation^ 
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Fig. 3. Effective index n e ff of the fundamental cladding 
mode versus normalized wavelength A/ A for holey fibers 
with normalized air-hole diameter d/A = 0.4 and varying 
base material. The dashed lines are the predictions of 
Eq. JSJl and the data points are results of fully-vectorial 
plane- wave simulations^ The inset shows the unit cell 
of the periodic cladding structure and the fundamental 
cladding eigenfunction ip c \ with 7^ = 15.5728, obtained 
with the aid of a finite-clement simulation. 7 



As expected the mode-field diameter decreases as the 
normalized air-hole diameter is increased and the mode 
becomes more localized. For Vpcf = 7T we find that 
io eff /A ~ 0.627. 

With Eqs. © and (JSJ at hand we have now provided 
a unified theory of the dispersion relation in the short- 
wavelength regime for PCFs with arbitrary base materi- 
als and Eq. © illustrates how geometrical confinement 
modifies the linear free-space dispersion relation. 

In fiber optics it is common to express the dispersion 
properties in terms of the effective index n G fr = c[3 /uj ver- 
sus the free-space wavelength A = c2it/lo. From Eq. © 
it follows straightforwardly that 



n c s = n b \ 



Y 



47T 2 7 



(8) 



which obviously is in qualitative agreement with the ac- 
cepted view that n e g increases monotonously with de- 
creasing wavelength and approaches n b in the asymptotic 
short-wavelength limit as reported for e.g. silica-based 
PCFs. 2 However, how good is the quantitative agreement 
for different base materials? In Figs. [21 and |21 we em- 



ploy fully-vectorial plane- wave simulations** to compare 
Eq. Q with the predictions of Eq. JHJ for the fundamen- 
tal core and cladding modes, respectively. For the core- 
modes we have employed a sufficiently large super-cell 
configuration. As seen there is an over-all good agree- 
ment between the fully-vectorial numerical results from 
Eq. JQ| and the semi- analytical predictions of Eq. In 
the short-wavelength limit A < A the agreement is excel- 
lent, which underlines the high relevance of the present 
results to large-mode area PCFs. As the index of the 
base material is increased the agreement is even better 
(results not shown). As the wavelength A is increased 
and becomes comparable to the pitch A the quantitative 
agreement is less good. The reason is well-known from 
small-core silica-based PCFs where a non-negligible frac- 
tion of the electrical field is forced to reside in the air-hole 
regions and where also vectorial effects of Eq. start 
to matter. 

In conclusion we have shown how a unified description 
of the short-wavelength dispersion and modal properties 
is possible. The theory illustrates how the waveguide dis- 
persion originates from the geometrical transverse local- 
ization of the mode and the semi-analytical description 
of the short-wavelength properties is readily applied to 
PCFs made from any base material. 
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